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Abstract 


Associations  between  cetacean  distributions,  oceanographic  features,  and  bioacoustic  backscatter  were  examined 
during  two  process  cruises  in  the  northern  California  Current  System  (CCS)  during  late  spring  and  summer  2000.  Line- 
transect  surveys  of  cetaceans  were  conducted  across  the  shelf  and  slope,  out  to  150  km  offshore  from  Newport,  Oregon 
(44.6°N)  to  Crescent  City,  California  (41 .9°N),  in  conjunction  with  multidisciplinary  mesoscale  and  fine-scale  surveys  of 
ocean  and  ecosystem  structure.  Occurrence  patterns  (presence/absence)  of  cetaceans  were  compared  with  hydrographic 
and  ecological  variables  (e.g.,  sea  surface  salinity,  sea  surface  temperature,  thermocline  depth,  halocline  depth, 
chlorophyll  maximum,  distance  to  the  center  of  the  equatorward  jet,  distance  to  the  shoreward  edge  of  the  upwelling 
front,  and  acoustic  backscatter  at  38,  120,  200  and  420  kHz)  derived  from  a  towed,  undulating  array  and  a  bioacoustic 
system.  Using  a  multiple  logistic  regression  model,  60.2%  and  94.4%  of  the  variation  in  occurrence  patterns  of 
humpback  whales  Megaptera  novaeangliae  during  late  spring  and  summer,  respectively,  were  explained.  Sea  surface 
temperature,  depth,  and  distance  to  the  alongshore  upwelling  front  were  the  most  important  environmental  variables 
during  June,  when  humpbacks  occurred  over  the  slope  (200-2000  m).  During  August,  when  humpbacks  concentrated 
over  a  submarine  bank  (Heceta  Bank)  and  off  Cape  Blanco,  sea  surface  salinity  was  the  most  important  variable, 
followed  by  latitude  and  depth.  Humpbacks  did  not  occur  in  the  lowest  salinity  water  of  the  Columbia  River  plume 
For  harbor  porpoise  Phocoena  phocoena,  the  model  explained  79.2%  and  70.1%  of  the  variation  in  their  occurrence 
patterns  during  June  and  August,  respectively.  During  spring,  latitude,  sea  surface  salinity,  and  thermocline  gradient 
were  the  most  important  predictors.  During  summer,  latitude  and  distance  to  the  inshore  edge  of  the  upwelling  front 
were  the  most  important  variables.  Typically  a  coastal  species,  harbor  porpoises  extended  their  distribution  farther 
offshore  at  Heceta  Bank  and  at  Cape  Blanco,  where  they  were  associated  with  the  higher  chlorophyll  concentrations  in 
these  regions.  Pacific  white-sided  dolphin  Lagenorhynchus  obliquidens  was  the  most  numerous  small  cetacean  in  early 
June,  but  was  rare  during  August.  The  model  explained  44.5%  of  the  variation  in  their  occurrence  pattern,  which  was 
best  described  by  distance  to  the  upwelling  front  and  acoustic  backscatter  at  38  kHz.  The  model  of  the' occurrence 
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pattern  of  Dali’s  porpoise  Phocoenoides  dalli  was  more  successful  when  mesoscale  variability  in  the  CCS  was  higher 
during  summer.  Thus,  the  responses  of  cetaceans  to  biophysical  features  and  upwelling  processes  in  the  northern  CCS 
were  both  seasonally  and  spatially  specific.  Heceta  Bank  and  associated  flow-topography  interactions  were  very 
important  to  a  cascade  of  trophic  dynamics  that  ultimately  influenced  the  distribution  of  foraging  cetaceans.  The  higher 
productivity  associated  with  upwelling  near  Cape  Blanco  also  had  a  strong  influence  on  the  distribution  of  cetaceans. 
©  2004  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  goal  of  US  GLOBEC  is  to  understand  and 
predict  how  marine  species  respond  to  global 
climate  change.  Among  the  uncertainties  in  a 
warmer  global  climate  is  the  extent  to  which 
upwelling  will  increase  or  decrease  in  specific 
boundary  current  systems,  such  as  the  California 
Current  System  (CCS),  and  consequently  affect 
the  productivity  and  structure  of  marine  commu¬ 
nities.  On  decadal  scales,  climate-related  shifts  in 
ecosystem  structure  of  the  CCS  have  been  well 
documented  (Roemmich  and  McGowan,  1995; 
Mantua  et  al.,  1997;  Peterson  and  Schwing,  2003). 
The  goal  of  our  research,  a  component  of  the  US 
GLOBEC  Northeast  Pacific  CCS  program  (Batch- 
elder  et  al.,  2002;  Strub  et  al.,  2002),  is  to 
determine  how  top-trophic  predators  (i.e.,  mam¬ 
mals  and  seabirds)  in  the  northern  CCS  relate  to 
middle  trophic  levels  and  bio-physical  coupling  in 
the  system.  Stemming  from  these  studies,  our  long¬ 
term  objective  is  to  develop  predictive  biophysical 
models  of  cetacean  occurrence  patterns  to  improve 
our  understanding  of  the  responses  of  top  pre¬ 
dators  to  climate-related  variability  in  the  struc¬ 
ture  of  an  upwelling  boundary  current  system. 

Enhanced  productivity  associated  with  coastal 
upwelling  systems  can  provide  important  predict¬ 
able  summer  foraging  for  large  whales  (Fiedler  et 
al.,  1998;  Gill,  2002).  The  frequency,  duration  and 
seasonal  cycle  of  upwelling  in  the  northern  CCS 
also  is  expected  to  have  important  effects  on  the 
productivity  and  community  structure  of  lower 
and  middle  trophic  levels,  and  hence  top-trophic 
levels.  During  strong  upwelling  years,  euphausiids 
are  the  dominant  food  consumed  by  many  species 
of  pelagic  nekton  in  the  CCS  (Ainley  and 
Boekelheide,  1990;  Brodeur  and  Pearcy,  1992; 
Ainley  et  al.,  1996).  In  central  California,  Ainley  et 
al.  (1996)  found  that  upwelling  conditions  are 


favorable  to  the  availability  of  the  euphausiid 
Thysanoessa  spinifera ,  an  important  prey  species 
for  many  fish,  seabirds  and  whales.  Off  southern 
California,  the  distribution  of  Balaenoptera  whales 
is  determined  by  their  attraction  to  areas  of 
predictable  high  densities  of  Euphausia  pacifica 
and  T.  spinifera  (Croll  et  al.,  1998).  During 
summer  and  fall,  blue  whales  Balaenoptera  muscu- 
lus  are  found  in  cold,  well-mixed,  productive  water 
that  upwells  along  the  coast  north  of  Point 
Conception  and  advects  south  (Fiedler  et  al., 
1998).  However,  in  more  northern  sections  of  the 
west  coast  of  North  America,  the  influence  of 
upwelling  and  dynamics  of  the  northern  CCS  on 
cetaceans  have  been  far  less  studied.  Our  GLO¬ 
BEC  study  is  the  first  multidisciplinary  program  in 
the  northern  CCS  to  apply  synoptic  integrative 
oceanographic  sampling  of  mesoscale  and  fine- 
scale  processes  to  the  study  of  cetacean  ecology. 

The  GLOBEC  Northern  CCS  study  region 
(41 .9-44. 6°N)  occurs  within  an  eastern  boundary 
current  system  that  extends  along  the  west  coast  of 
North  America  from  the  Strait  of  Juan  de  Fuca  to 
the  tip  of  Baja  California  (Hickey,  1998).  Nor¬ 
mally,  during  the  spring  and  summer  upwelling 
season,  mean  wind  stresses  are  southward,  sea 
levels  are  low,  and  sea  surface  temperature  (SST)  is 
cool  over  the  northeast  Pacific  continental  shelf 
(Strub  et  al.,  1987).  When  the  upwelling  index  is 
high,  an  alongshore,  subsurface  ribbon  of  cool 
water  and  an  equatorward  surface  jet  are  usually 
observed  over  the  shelf  (Huyer  and  Smith,  1974; 
Smith  et  al.,  1999).  On  a  finer  scale,  during 
persistent  upwelling,  broad  tongues  and  narrow 
filaments  of  cold,  chlorophyll-rich  water  may 
extend  over  100  km  offshore  and  downstream, 
usually  associated  with  coastal  promontories 
(Hood  et  al.,  1990;  Barth  et  al.,  2002).  The  latter 
applies  to  Cape  Blanco,  an  important  topographic 
feature  in  our  study  region.  Here  the  upwelling  jet 
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separates  from  the  shelf  and  meanders  equator- 
ward  as  an  oceanic  jet  (Strub  et  al.,  1991;  Batteen, 
1997;  Barth  et  al.,  2000).  To  examine  the  ecological 
importance  of  the  flow-topography  interactions  at 
Cape  Blanco  on  productivity,  and  consequently  on 
cetacean  distribution,  is  among  our  objectives.  In 
addition,  our  GLOBEC  study  was  designed  to 
examine  the  ecological  influence  of  one  of  the 
major  banks  in  the  CCS  (Heceta  Bank, 
44^I4.5°N).  We  hypothesize  that  upwelling,  flow- 
topography  interactions,  and  local  recirculations 
associated  with  the  bank  generate  enhanced 
seasonal  productivity  in  this  region,  and  that  this 
will  attract  foraging  cetaceans. 


2.  Materials  and  methods 

2. 1.  Oceanographic  data 

Hydrographic  data  were  collected  using  a 
towed,  undulating  vehicle  known  as  ‘SeaSoar’ 
(Pollard,  1986),  cycled  rapidly  from  the  surface  to 
depth  while  being  towed  at  3.6 ms-1  (7 kts).  The 
vehicle  was  equipped  with  a  Seabird  SBE  911  + 
conductivity-temperature-depth  (CTD)  instrument 
with  pumped,  dual  temperature-conductivity  sen¬ 
sors  pointing  forward  through  a  hole  in  the 
SeaSoar  nose.  A  Western  Environmental  Technol¬ 
ogy  Laboratories  (WET  Labs)  Flashpak  fluorom- 
eter,  using  green  excitation  (490,  30  nm  bandpass) 
and  685-nm  detection  wavelength,  was  used  to 
estimate  chlorophyll  fluorescence.  The  signal  was 
converted  to  chlorophyll  (mg  m-3)  using  a  calibra¬ 
tion  developed  from  discrete  samples  measured 
with  high-performance  liquid  chromatography. 
Most  CTD-data  were  collected  by  cycling  SeaSoar 
on  a  bare  hydrographic  cable  from  the  sea  surface 
to  1 15- 120  m  over  deep  water  and  to  within  10  m 
of  the  bottom  over  the  shelf.  Cycle  time  in  deep 
water  was  about  6.5  min.  The  result  was  hydro- 
graphic  data  with  high  spatial  resolution  (1.3  km 
between  surface  points  and  500  m  between  profiles 
at  mid-depth)  obtained  rapidly  (cross-margin 
sections  in  2-10  h  and  large-area  maps  in  2-6 
days)  so  that  a  detailed,  synoptic  view  of  the 
system  was  available.  Over  the  shallow  continental 
shelf,  SeaSoar  was  cycled  from  0-55-0  m  every 


1.5  min,  which  resulted  in  profiles  separated  by 
about  300  m  at  the  surface  and  170  m  at  mid¬ 
depth.  In  order  to  sample  to  greater  depths  over 
the  continental  slope,  SeaSoar  also  was  towed 
using  a  faired  cable  cycling  from  0-300-0  m  every 
8  min.  Concurrently,  data  were  logged  from  a  hull- 
mounted  Acoustic  Doppler  Current  Profiler 
(ADCP).  Further  details  of  this  hydrographic  data 
set  can  be  found  in  Barth  et  al.  (2005). 

Using  time  varying  lags  and  an  optimized 
thermal  mass  correction,  the  24-Hz  temperature 
and  conductivity  data  were  realigned  and  cor¬ 
rected  to  calculate  24-Hz  salinity,  and  averaged  to 
yield  1-Hz  values  (Barth  et  al.,  2000).  The  final 
1-Hz  data  files  contain  unfiltered  GPS  latitude 
and  longitude;  pressure;  temperature;  salinity; 
and  density  anomaly  (sigma-/,  for  example  as 
1024  kgm-3  minus  1000)  computed  using  the  1980 
equation  of  state;  chlorophyll  concentration 
(mgm-3);  date  and  time.  The  1-Hz  SeaSoar  CTD 
data  were  first  averaged  vertically  into  2-dbar  bins, 
from  which  environmental  indices  were  derived. 
These  included  sea  surface  temperature  (SST 
where  surface  refers  to  5-m  depth),  sea  surface 
salinity  (SSS),  thermocline  depth,  halocline  depth, 
and  chlorophyll  maximum.  The  indices  were 
averaged  into  along-track  bins  of  approximately 
15  min  (or  4.4  km  at  average  ship  speed  of  9.5 
knots)  for  comparison  with  seabird  (Ainley  et  al., 
2005)  and  mammal  data. 

Acoustic  volume  backscattering  (dB)  data  were 
collected  with  a  Hydroacoustics  Technology, 
Incorporated  (HTI)  model  244  echosounder  using 
frequencies  of  38,  120,  200,  and  420kHz.  The  four 
transducers  were  deployed  on  a  towed  body  at 
about  4  m  depth  off  the  port  side  of  the  ship,  in 
combination  with  the  SeaSoar  vehicle  off  the  stern. 
The  HTI  instrument  collected  echo  integration 
data  in  1  m  bins  using  a  raw  ping  rate  of  4 pings-1 
or  1  pings-1  for  each  frequency.  These  data  were 
averaged  into  12-s  ensembles,  resulting  in  a 
horizontal  distance  of  about  50  m  between 
adjacent  ensembles.  For  this  study,  we  use  depth- 
integrated  backscatter  from  the  surface  to  100  m, 
and  averaged  over  15-min  segments  of  trackline, 
for  comparison  with  seabird  and  cetacean 
data.  Roughly,  the  38  kHz  backscatter  gives  an 
approximate  measure  of  fish,  the  120  and  200  kHz 
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can  be  associated  with  larger  and  smaller  euphau- 
siids  respectively,  and  the  420  kHz  is  most 
closely  related  to  copepods.  Further  discussion 
of  the  interpretation  of  bioacoustic  data  can  be 
found  in  Ressler  et  al.  (2005),  and  references 
therein. 

The  spatially  averaged  temperature,  salinity  and 
pressure  data  were  used  to  compute  geopotential 
anomaly  (dynamic  height  in  meters  multiplied  by 
the  acceleration  of  gravity)  in  Jkg-1  (m2s~2) 
relative  to  lOOdb.  On  east-west  sections  where 
SeaSoar  profiles  were  shallower  than  lOOdb, 
geopotential  anomaly  was  calculated  using  the 
extrapolation  technique  described  by  Reid  and 
Mantyla  (1976).  The  geopotential  anomaly  of 
2.35  m2s-2  is  indicative  of  the  center  of  the 
equatorward  jet  and  the  2.0  nr  s  2  marks  the 
inshore  side  of  the  jet.  Distances  to  each  of  these 
features  were  included  as  variables  in  the  analyses. 
In  addition,  since  the  dynamic  height  reflects  the 
subsurface  structure,  we  also  tested  distance  to  a 
surface  feature  (i.e.  distance  to  the  SST  isotherm 
11.5°C  during  June  and  12.0°C  during  August) 
that  best  represented  the  inshore  side  of  the 
alongshore  upwelling  front. 


2.2.  Cetacean  surveys 

Line-transect  surveys  of  cetaceans  were  con¬ 
ducted  during  the  periods  May  29-June  13  and 
July  27-August  12,  2000  across  the  shelf  and 
slope  off  Oregon  and  northern  California 
(41 .9-44. 65°N)  from  the  coast  to  ~  1 25. 5°W 
(Fig.  1).  Surveys  were  conducted  in  passing  mode 
while  the  R.V.  New  Horizon  was  in  transit  between 
stations  for  hydrographic  and  zooplankton  sam¬ 
pling.  Surveys  were  designed  to  follow  12  mesos- 
cale  across-shelf  and  slope  transects  and  additional 
fine-scale  survey  lines  similar  to  those  run  by  a 
second  ship,  R.V.  Wecoma,  for  acquisition  of 
hydrographic  data  (see  Barth  et  al.,  2005;  their 
Fig.  1).  Observations  were  conducted  from  the 
flying  bridge  during  daylight  (M)600-2030)  when¬ 
ever  sufficient  visibility  (i.e.,  ^2nm)  and  weather 
(i.e.,  ^Beaufort  6)  allowed.  Two  observers  simul¬ 
taneously  surveyed  to  the  horizon  with  25  x  1 50 
binoculars,  equipped  with  compass  and  reticle. 
Each  observer  surveyed  a  1 00°  arc  from  1 0  off  the 
bow  (opposite  side)  to  90°  on  their  side  of  the  ship. 
A  third  observer  focused  on  the  track-line  by  eye, 
aided  with  7  x  50  hand-held  binoculars.  Sightings 


126°  125°  124' 


126°  125°  124° 


Fig.  1.  Cetacean  survey  effort  using  25  x  150  binoculars  (solid  bold  line  segments)  during  May  29-June  13  and  during  July  27  August 
12,  2000  in  the  northern  California  Current  System. 
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were  entered  immediately  on  a  laptop  computer 
connected  to  the  ship's  GPS  system.  Positions  of 
all  sightings  were  corrected  to  reflect  the  actual 
location  of  the  cetaceans,  rather  than  the  ship's 
position.  The  height  from  the  surface  of  the  water 
to  the  eyes  of  observers  on  the  25  x  150  binoculars 
was  10.87  m.  Under  optimal  conditions,  balaenop- 
terids  observed  up  to  8.0  km  from  the  ship  could 
be  identified  to  species;  whales  on  the  horizon 
(10.4  km)  could  not  be  identified  to  species. 
Abundance  estimation  and  related  methodology 
will  be  reported  elsewhere. 

2.3.  Modeling  occurrence  pal  terns  of  cetaceans 

Occurrence  (presence/absence)  patterns  of  ceta¬ 
ceans  relative  to  ocean  processes  and  biological 
backscatter  were  examined  for  four  species  having 
a  sufficiently  large  sample  size  (i.e.  number  of 
sightings  obtained  while  on  25  x150  binocular 
effort):  humpback  whale  Megaptera  novaeangliae , 
harbor  porpoise  Phocoena  phocoena ,  Dali  s  por¬ 
poise  Phocoenoides  dalli ,  and  Pacific  white-sided 
dolphin  Lagenorhynchus  obliquidens.  Each  sighting 
for  these  four  species  was  assigned  to  the  nearest 
15-min  segment  of  track-line  (or  approximately 
4.4  km  at  average  ship  speed  of  9.5  kts)  for  which  a 
suite  of  16  hydrographic,  acoustic  and  ecological 
variables  were  available  (see  above).  The  time 
difference  between  the  collection  of  hydrographic 
data  and  the  acquisition  of  line-transect  survey 
data  was  never  more  than  about  1 2  h,  as  the  two 
ships  were  never  separated  by  more  than  this  time 
lag.  The  coherence  in  time  and  space  of  the  inshore 
(depth  <-100m)  and  offshore  (depth  >-  100  m) 
circulation  structures  are  typically  greater  than  the 
distances  between  cetaceans  and  the  track-line 
(Barth  et  al.,  2005).  For  example,  the  offshore 
meanders  are  long-lived  and  last  for  weeks  to 
months.  Inshore,  currents  throughout  the  water 
column  and  surface  temperature  and  salinity 
change  with  the  wind  forcing  on  a  time  scale 
2  10  days.  The  deeper  (^50  m)  horizontal  density, 
temperature  and  salinity  fronts  are  more  stable.  In 
general,  the  alongshore  correlation  length  scales 
are  much  longer  than  the  cross-shelf  correlation 
length  scales.  Kundu  and  Allen  (1976),  using 
moored  velocity  data,  found  alongshore  and 


cross-shelf  correlation  scales  of  at  least  30  km 
(i.e.  about  the  same  size  as  our  distances  between 
mesoscale  survey  lines,  and  significantly  larger 
than  the  distances  of  cetacean  sightings  to  the 
track-line  and  the  nearest  grid  segment  of  SeaSoar- 
derived  environmental  data). 

Multiple  logistic  regression  was  used  to  model 
cetacean  occurrence  (Cox  and  Snell,  1989),  be¬ 
cause  outcomes  were  treated  as  binary  (a  cetacean 
species  was  either  present  or  absent  within  a  15- 
min  transect).  An  important  assumption  of  logistic 
regression  is  that  outcomes  are  binomially  dis¬ 
tributed  (Cox  and  Snell,  1989);  therefore,  we 
confirmed  that  occurrence  of  each  species  did  not 
deviate  from  binomial  distribution  (goodness  of  fit 
X~-test,  P>0.05).  We  also  tested  for  autocorrela¬ 
tion  of  cetacean  occurrence  among  the  15-min 
transect  segments  for  each  species  using  (the  ‘acf 
command  of)  S-PLUS  (S-PLUS,  1997).  Occur¬ 
rence  patterns  were  not  autocorrelated  among 
transects  for  any  of  the  species. 

Maximum-likelihood  logistic  regression  was 
implemented  using  program  STATA  (STATA 
Corporation,  1995);  the  logit'  command  for  the 
likelihood  ratio  test  was  used  to  examine  the 
occurrence  of  each  species  by  season  (i.e.  May 
June  and  July  August)  relative  to  two  spatially 
fixed  variables  (latitude,  depth)  and  the  following 
14  spatially  varying  variables:  SST,  SSS,  depth  of 
the  thermocline,  thermocline  gradient,  depth  of  the 
halocline,  halocline  gradient,  distance  to  the  jet 
‘center'  (geopotential  anomaly  of  2.35  m2s-2), 
distance  to  the  inshore  edge  of  the  jet  (geopotential 
anomaly  of  2.0  irrs  “),  distance  to  the  inshore 
edge  of  the  surface  front  (11.5°C  SST  isotherm 
during  spring  and  12.0  C  SST  isotherm  during 
summer),  value  of  chlorophyll  a  maximum,  and 
integrated  acoustic  backscatter  at  four  frequencies 
(38,  120,  200,  and  420kHz).  A  total  of  420  grid- 
segments,  for  which  cetacean  effort  (while  on 
25  x  150  binoculars)  and  the  full  suite  of  SeaSoar- 
derived  and  bioacoustic  data  were  available,  were 
used  in  the  multiple  logistic  regression  models:  226 
segments  in  May-June  and  194  segments  in 
July-August  (Fig.  2).  We  report  the  likelihood 
ratio  statistic,  which  is  equivalent  to  ‘deviance'  as 
reported  in  generalized  linear  models,  and  there¬ 
fore  is  analogous  to  sums  of  squares  reported  in 
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Fig.  2.  Location  of  the  starting  position  for  each  15-min  survey  segment  (square)  used  in  the  multiple  logistic  regression  models  of 
cetacean  occurrence,  May-June  and  July-August  2000,  relative  to  chlorophyll  a  concentration  at  5  m.  All  segments  used  in  the  models 
had  coverage  of  line-transect  survey  effort  for  cetaceans,  oceanographic  sampling  using  the  towed  SeaSoar  vehicle,  and  bioacoustic 
sampling  using  the  multi-frequency  system. 


ANOVA.  We  evaluated  departure  from  linearity 
by  testing  whether  the  coefficients  of  a  polynomial 
logistic  regression  model  (to  order  2)  differed 
significantly  from  zero. 

We  modeled  occurrence  patterns  for  each 
species  by  season,  except  for  Pacific  white-sided 
dolphin.  The  number  of  sightings  (3)  of  Pacific 
white-sided  dolphins  during  August  2000  was 
insufficient  for  meaningful  analysis;  therefore,  we 
modeled  occurrence  pattern  for  this  species 
after  combining  spring  and  summer  data.  To 
control  for  possible  diurnal  patterns  in  cetacean 
behavior  or  migratory  activity  during  the  cruise 
period,  we  included  time  of  day  and  Julian  date 
in  the  analyses.  In  addition,  wind  speed  was 


included  in  an  initial  model  to  check  for  possible 
effects  of  occurrence  pattern  related  to  either 
animal  activity  during  increased  sea  state  or 
reduced  probability  of  detection  with  increased 
sea  state.  Wind  speed  was  insignificant  for  5  of  the 
7  models  (humpback  whales  in  both  seasons, 
Dali's  porpoise  in  both  seasons,  and  harbor 
porpoise  in  spring)  and  adding  wind  speed  reduced 
the  explained  variance  in  every  model  (except  for  a 
slight  improvement  in  the  model  of  Pacific  white¬ 
sided  dolphin).  Therefore,  to  develop  the  best 
predictive  models  and  to  test  the  direct  effects  of 
ocean  conditions  resulting  from  wind-driven  up- 
welling,  wind  speed  was  not  included  in  the  final 
models. 
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The  multiple  logistic  regression  analyses  were 
conducted  by  initially  entering  all  16  independent 
variables  together  in  each  model.  Using  a  stepwise 
procedure,  insignificant  terms  were  dropped,  one 
at  a  time,  in  order  of  increasing  P-values.  Because 
many  terms  were  correlated  (Table  2),  importance 
of  some  variables  was  likely  masked  by  others  in 
the  initial  model.  Therefore,  we  tested  for  effects  of 
eliminated  terms  by  putting  them,  one  at  a  time, 
back  in  the  model.  The  model  was  complete  when 
no  terms  could  be  added  or  dropped.  In  a  multiple 
regression  model,  any  independent  variables  that 
test  as  having  a  significant  relationship  with  the 
dependent  variable  (cetacean  occurrence)  are  true 
influences  (Seber,  1977;  Kleinbaum  et  al.,  1988). 
That  is,  their  effects  are  independent  of  those  of 
other  independent  variables  (correlated  or  not) 
that  also  have  significant  relationships  with 
occurrence,  because  each  independent  term  in¬ 
cluded  in  the  model  is  evaluated  while  taking  into 
account  (controlling  for)  the  effect  of  each  of  the 


other  independent  variables  (see  Ainley  et  al., 
2005,  for  additional  detail  on  the  properties  and 
benefits  of  multiple  regression  analyses  when 
applied  to  correlated  variables). 


3.  Results 

3.1.  Correlations  among  physical  and  biological 
oceanographic  variables 

The  physical  and  biological  environmental 
variables  were  correlated  with  one  another  in  99 
(73%)  of  the  136  possible  correlations  (Table  1). 
Considering  the  relationships  having  r-values 
>0.4,  (P< 0.0001),  for  reasons  of  brevity,  the 
following  was  indicated:  (1)  distance  to  the  center 
of  the  equatorward  jet,  sea  surface  salinity,  and 
abundance  of  smaller  prey  (represented  by 
>  120  kHz  backscatter)  increased  with  Julian  date, 
while  density  of  larger  prey  (38  kHz  size  range) 


Table  1 

The  Pearson  linear  correlation  coefficients  (r)  among  physical  and  biological  environmental  variables51;  n  =  420  survey  transects 
up  welling  period  2000  J 


JP  LT  FRT  JTI  JTC  DPT  SST  SSS  TDP  TST  HDP  HST  CM  38kHz  120kHz  200kHz 


LT 

0.03 

FRT 

0.04 

-0.11 

JTI 

0.18 

0.04 

0.80 

JTC 

0.57 

-0.03 

0.33 

0.41 

DPT 

0.23 

-0.31 

0.66 

0.72 

0.21 

SST 

-0.07 

-0.25 

0.49 

0.44 

-0.21 

0.58 

SSS 

0.70 

-0.23 

-0.09 

-0.04 

0.57 

-0.01 

-0.43 

T-DP 

-0.29 

-0.24 

0.40 

0.33 

-0.15 

0.31 

0.43 

-0.23 

T-ST 

0.04 

-0.06 

0.29 

0.25 

-0.23 

0.36 

0.78 

-0.34 

0.20 

H-DP 

-0.05 

-0.28 

0.33 

0.27 

0.06 

-0.34 

0.37 

-0.02 

0.61 

0.20 

H-ST 

-0.71 

0.05 

-0.06 

-0.14 

-0.49 

-0.13 

0.24 

-0.77 

0.10 

0.26 

0.01 

CM 

0.33 

0.07 

-0.21 

-0.17 

0.38 

-0.31 

-0.32 

0.40 

-0.30 

-0.19 

-0.23 

-0.25 

38  kHz 

-0.55 

-0.22 

0.20 

0.05 

-0.32 

0.17 

0.28 

-0.43 

0.25 

0.18 

0.11 

0.38 

-0.25 

120kHz 

0.40 

0.00 

0.21 

0.25 

0.19 

0.37 

0.20 

0.14 

-0.07 

0.21 

0.00 

-0.28 

0.07 

0.38 

200  kHz 

0.57 

0.20 

0.14 

0.33 

0.34 

0.35 

0.05 

0.25 

-0.15 

0.16 

-0.04 

-0.37 

0.13 

-0.03 

420  kHz 

0.47 

0.08 

0.30 

0.37 

0.34 

0.43 

0.15 

0.14 

-0.15 

0.16 

-0.05 

-0.27 

0.00 

-0.02 

0.72 


■'Relationships  having  r-values  in  bold  were  significantly  correlated  (/><0.05).  LT,  latitude;  JD,  Julian  date;  FRT.  distance  to 

upwelling  surface  front  defined  by  1 1 .5  C  surface  isotherm  during  June  and  12.0  C  surface  isotherm  during  August;  JTI  distance  to 
the  inshore  edge  of  the  jet  (geopotential  anomaly  of  2.0  mV2);  JTC,  distance  to  the  jet  ‘center’  (geopotential  anomaly  of  2  35  mV2)* 
DPT  ocean  depth;  SST,  sea  surface  temperature;  SSS.  sea  surface  salinity;  TDP,  thermocline  depth;  TST,  thermocline  strength;  HDp[ 
ha  Iodine  depth;  HST,  halochne  strength;  CM,  value  of  chlorophyll  maximum;  38  kHz,  integrated  backscatter  O-IOOm;  120  kHz 
backscatter;  200kHz  backscatter;  and  420kHz  backscatter.  A  positive  correlation  between  acoustic  backscatter  and  another  variable 
indicates  a  positive  relationship,  and  the  reverse  for  a  negative  correlation. 
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decreased  with  date:  (2)  SST  and  ocean  depth 
increased  with  increase  in  distance  between  SST- 
defined  shoreward  edge  of  upwelling  front  and 
dynamic  height  defined  shoreward  edge  of  jet  (the 
latter  also  were  highly  correlated);  (3)  SSS  increased 
and  halocline  strength  decreased  with  increase  in 
distance  from  the  center  of  the  equatorward  jet;  (4) 
SST  and  density  of  the  smallest  prey  size  class 
increased  with  increase  in  ocean  depth;  (5)  thermo- 
cline  depth  and  gradient  increased,  while  SSS 
decreased,  with  increase  in  SST;  (6)  value  of  the 
chlorophyll  maximum  increased,  halocline  strength 
and  abundance  of  larger  prey  (38  kHz)  decreased 
with  increase  in  SSS;  (7)  halocline  depth  increased 
with  increase  in  thermocline  depth;  (8)  backscatter 
at  the  three  highest  frequencies  increased  in  accord 
with  one  another. 

3.2.  Cetacean  species  composition 

During  August,  the  cetacean  survey  coverage 
was  less  uniform  across  the  shelf  at  43.5  N  than 
during  June,  in  order  that  the  biophysical  struc¬ 
ture  on  the  outside  of  a  warm  mesoscale  eddy 
also  could  be  studied  (Fig.  1).  Survey  coverage 


extended  west  over  the  basin  to  126  W  to  reach  the 
outside  of  this  feature.  A  total  of  223  sightings  of 
cetaceans  (1698  animals)  were  recorded  during 
spring  and  298  sightings  (614  animals)  were 
recorded  during  summer  (Table  2).  Sixteen  ceta¬ 
cean  species  were  recorded  in  the  northern  CCS 
during  2000  (Table  2):  four  species  of  large  whale 
(humpback  M.  novaeangliae ,  gray  Eschrichtius 
robustus ,  fin  Balaenoptera  physalus ,  and  sperm 
whale  Physeter  macrocephalus );  seven  species  of 
medium-sized  whale  (minke  Balaenoptera  acutor - 
ostrata ,  killer  whale  Oreinus  orca ,  Baird's  beaked 
whale  Berardius  bairdii ,  Cuvier's  beaked  whale 
Ziphius  cavirostris ,  Mesoplodon  sp..  Hyper oodon 
sp.,  and  Risso's  dolphin  Grampus  griseus ,),  and 
five  species  of  small  cetacean  (Pacific  white-sided 
dolphin  L.  obliquidens ,  northern  right  whale 
dolphin  Lissodelphis  borealis .  Dali's  porpoise  P. 
dalli ,  harbor  porpoise  Phocoena  phoeoena ,  and 
dwarf  sperm  whale  Kogia  simus). 

3.3.  Occurrence  patterns  of  cetaceans 

The  multiple  logistic  regression  models 
explained  25.6  94.4%  of  the  variance  in  the 


Table  2 

Summary  of  cetacean  sightings  and  numbers  of  animals  observed  during  line-transect  surveys  conducted  in  conjunction  with  two 
GLOBEC  Northeast  Pacific  Northern  California  Current  cruises  during  2000 


Species 

Survey  1:  May  29-June  13 

Survey  2:  July  27  August  13 

#  Sightings 

#  Animals 

#  Sightings 

#  Animals 

Pacific  white-sided  dolphin 

86 

1063 

3 

40 

Dali’s  porpoise 

31 

99 

130 

294 

Humpback  whale 

35 

70 

45 

76 

Harbor  porpoise 

41 

62 

87 

125 

Northern  right  whale  dolphin 

9 

195 

0 

0 

Risso’s  dolphin 

7 

87 

3 

17 

Killer  whale 

4 

19 

2 

9 

Sperm  whale 

2 

6 

3 

3 

Baird's  beaked  whale 

5 

27 

3 

13 

Mesoplodon  sp. 

0 

0 

4 

6 

Cuvier’s  beaked  whale 

0 

0 

1 

1 

Hyperoodon  sp. 

0 

0 

1 

3 

Gray  whale 

1 

1 

6 

14 

Fin  whale 

0 

0 

6 

9 

Minke  whale 

1 

1 

4 

4 

K.  simus 

1 

1 

0 

0 
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Table  3 


Results  of  multiple  logistic  regression  models  to  examine  occurrence  patterns  of  Pacific  white-sided  dolphin  (PWSD).  humpback  whale 
(HBW)  Dali  s  porpoise  (DAPO)  and  harbor  porpoise  (HAPO)  relative  to  physical  and  biological  oceanographic  variables  during 
spring  (late  May-early  June)  and  summer  (late  July-early  August)  2000 


PWSD 

HBWspr 

HBWsum 

DAPOspr 

DAPOsum 

HAPOspr 

HAPOsum 

%  Variation  explained 

44.5 

60.2 

94.4 

25.6 

46.1 

79.2 

70.1 

Latitude 

_ 2 

+  i 

.  i 

Distance  to  Front  A 

— 

3 

_ 

M 

+ 

*r 
_ 2 

Distance  to  2.0  m2s-2 

M1 

+ 

+ 

_ 2 

Distance  to  2.35  m2s~2 

_ 

+  3 

Ocean  depth 

+ 

+  2 

3 

+  1 

_ 

SST 

+  * 

_ 

M 

+ 

SSS 

- 

M 

M1 

+  2 

Thermocline  Depth 

_ 

Thermocline  Gradient 

M3 

Halocline  depth 

Halocline  gradient 

M 

Backscatter  38  kHz 

A  +  2 

c  + 

B  + 

A  +  3 

120kHz 

B _ 3 

CM 

A  + 

200  kHz 

B  + 

M 

B _ 2 

CM 

B_ 

420  kHz 

A  +  « 

A  + 

CM 

Chlorophyll  maximum 

C_ 

A  + 

A  + 

B  + 

Shown  are  the  signs  of  the  regression  coefficients  for  variables  significantly  related  to  occurrence:  ( + ),  positive  response  to  increase  in  a 

variable;  (-),  negative  response;  and  “M  \  curvilinear  response  where  a  significant  linear  effect  was  not  observed  Superscripts  (1  2 
and  3)  adjacent  to  a  coefficient  indicate  that  the  variable  explained  the  first,  second,  or  third  largest  amount  of  variance  in  a  species*’ 
occurrence;  superscripts  (A,  B.  and  C)  indicate  that  the  variable  explained  the  first,  second,  or  third  largest  amount  of  variance  among 
the  five  biological  variables  (acoustic  backscatter  at  four  frequencies,  and  value  of  chlorophyll  a  maximum).  Details  of  the  variables  are 
presented  in  the  methods  section:  distance  to  the  jet  ‘center*  (geopotential  anomaly  of  2.35  m2s-2);  distance  to  the  inshore  edge  of  the 
jet  (geopotential  anomaly  of  2.0m  s--);  and  distance  to  inshore  surface  expression  of  upwelling  front  (1 1.5°C  SST  isotherm  during 
June  and  12.0  C  SST  isotherm  during  August).  The  latter  surface  front  is  referred  to  as  Front  A  above.  A  negative  coefficient  for 
distance  to  an  ocean  feature  indicates  that  the  species  occurred  near  the  feature.  A  positive  coefficient  for  cetacean  response  to 
backscatter  indicates  that  the  species  occurred  more  frequently  at  higher  densities  of  the  plankton  or  nekton  described  acousticallv  bv 
that  backscatter.  J 


distributions  of  four  cetacean  species  by  season 
(Table  3).  The  number  of  variables  significantly 
related  to  occurrence  pattern  ranged  from  two  (of 
the  possible  16  variables)  for  Dali's  porpoise 
(spring)  to  15  significant  variables  for  harbor 
porpoise  (summer).  The  variable  most  often 
significantly  related  to  occurrence  pattern  was 
the  distance  to  the  inshore  edge  of  the  upwelling 
front  (SST-defined).  This  variable  was  significant 
in  six  of  the  seven  models  (Table  3).  Other 
important  variables  that  were  significant  in  four 
or  five  models  included:  distance  to  the  inshore 
edge  of  the  jet,  distance  to  the  center  of  the  jet, 
ocean  depth,  sea  surface  salinity,  acoustic  back¬ 


scatter  at  three  frequencies  (38,  120  and  200  kHz), 
and  value  of  the  chlorophyll  maximum. 

3.3.1.  Humpback  whale 

The  humpback  whale  was  the  most  abundant 
large  cetacean  during  late  spring  and  summer 
(Table  2),  although,  their  distribution  varied 
greatly  with  season  (Fig.  3).  Spatially,  humpback 
whales  exhibited  a  bimodal  cross-shelf  distribution 
in  spring.  Some  humpbacks  occurred  in  cold, 
saline,  upwelled  waters  of  higher  surface  chlor¬ 
ophyll  north  and  south  of  Cape  Blanco  (Fig.  2). 
Most  humpbacks,  however,  remained  over  the 
slope  (200  2000  m)  in  warmer  SSTs  (>  12  °C)  and 
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intermediate  surface  salinities  (Fig.  3).  A  multiple 
logistic  regression  model  explained  60.2%  of  the 
variance  in  their  distribution  pattern  during  spring 
(Table  3).  Humpback  distribution  was  most 
correlated  with  sea  surface  temperature,  depth, 
and  distance  to  the  surface  expression  of  the 
upwelling  front  (i.e.  whales  were  near  the  11.5°C 
surface  isotherm).  In  addition,  humpback  occur¬ 
rence  was  significantly  negatively  correlated  with 
distance  to  the  kjet'  center  (i.e.  whales  over  the 
slope  are  near  the  geopotential  anomaly 
2.35  mss“2)  and  thermocline  depth  (i.e.  whales 
were  associated  with  shallower  mixed  layer)  (Table 
3).  Humpback  occurrence  was  significantly  posi¬ 
tively  correlated  with  distance  to  the  shoreward 
edge  of  the  jet  (i.e.  whales  were  farther  from  the 
geopotential  anomaly  2.0  mss-2),  integrated  back- 
scatter  at  38  kHz  (i.e.  larger  acoustic  targets  such 
as  fish)  and  200  kHz  (i.e.  zooplankton  such  as 
euphausiids),  as  well  as  the  value  of  the  chlor¬ 
ophyll  maximum  (Table  3).  Of  the  latter  three 
biological  variables,  value  of  the  chlorophyll 
maximum  was  most  strongly  correlated  with 
humpback  occurrence.  During  spring,  humpback 
occurrence  had  a  significant  curvilinear  (nonlinear) 
relationship  to  surface  salinity.  Humpback  whales 
did  not  occur  in  the  lowest  salinity  water  ( <  30.8) 
of  the  Columbia  River  plume  (Fig.  3). 

The  multiple  logistic  regression  model  of  hump¬ 
back  whale  distribution  during  August  explained 
94.4%  of  the  variance  in  occurrence  pattern  (Table 
3).  Surface  salinity  was  the  most  significant 
variable,  followed  by  latitude  and  depth.  As  was 
the  case  in  spring,  humpback  occurrence  during 
summer  had  a  significant  curvilinear  (nonlinear) 
relationship  to  surface  salinity,  and  humpbacks 
appear  to  avoid  the  low-salinity  surface  plume 
offshore  (Fig.  3).  The  negative  correlation  with 
both  distance  to  the  inshore  edge  of  the  surface 
front  (11.5°C  isotherm)  and  distance  to  the  jet 


center  reflect  the  bimodal  distribution  of  hump¬ 
backs  in  the  northern  region  (-44  N)  (Table  3  and 
Fig.  3).  Although  most  humpbacks  occurred  on 
the  shelf  near  the  shoreward  edge  of  the  upwelling 
front,  a  few  whales  occurred  on  the  slope 
associated  with  the  strong  mesoscale  variability 
in  the  jet.  The  seasonal  shift  in  humpback 
distribution  to  the  broad  shelf  at  Heceta  Bank 
(44  44.5  N)  is  reflected  in  the  significant  negative 
correlation  with  depth.  Humpback  whales  on  the 
shelf  typically  occurred  in  depths  <100m  on 
Heceta  Bank  and  off  Cape  Blanco  (Fig.  3).  Of 
the  biological  variables  tested  in  the  model,  value 
of  the  chlorophyll  maximum  was  most  strongly 
correlated  with  humpback  occurrence,  followed  by 
integrated  backscatter  at  38  kHz  (Table  3  and  Fig. 
4).  Humpback  whales  occurred  in  upwelled  waters 
on  Heceta  Bank  and  off  Cape  Blanco,  associated 
with  higher  surface  chlorophyll,  higher  integrated 
backscatter  at  38  kHz  and  high  densities  of 
juvenile  salmon  (Fig.  5).  On  Heceta  Bank,  hump¬ 
backs  occurred  within  a  zone  of  lower  mean  near¬ 
surface  velocities  and  recirculation  (see  ADCP- 
derived  velocity  vectors  at  5  m.  Fig.  6). 

3.3.2.  Other  balaenopterids 

Fin  whales  were  only  found  during  August, 
when  survey  effort  extended  further  offshore 
1 26  W;  >  2000  m)  to  the  western  edge  of  an 
anticyclonic  mesoscale  eddy  or  meander  (43.75°N) 
(Fig.  6).  Fin  whales  occurred  in  the  surface  flow  of 
higher  chlorophyll  that  was  advected  northward 
on  the  western  side  of  the  meander  (see  surface 
chlorophyll.  Fig.  4).  Fin  whales  appear  to  prefer 
the  outer,  basin  region  of  influence  of  the  coastal 
jet.  Therefore,  during  summer,  fin  whales  and 
humpback  whales  are  widely  separated  in  the 
northern  CCS.  Gray  whales  were  observed  very 
close  to  the  coast  in  cold  (<9.5°C)  upwelled 
waters,  especially  near  Cape  Blanco  in  August 


◄ — - - 

Fig.  3.  The  correspondence  of  humpback  whale  sightings  (*)  with  sea  surface  (5m)  temperature,  surface  salinity,  and  geopotential 
anomaly  at  5  m  relative  to  100  m,  during  May-June  and  July-August,  2000.  Several  features  related  to  upwelling  were  included  in  the 
models  of  cetacean  occurrence:  the  center  of  the  equatorward  jet,  defined  by  geopotential  anomaly  2.35  mV2  (the  double  heavy  line)- 
the  shoreward  edge  of  the  jet,  defined  by  geopotential  anomaly  2.0  mV2  (first  single  heavy  line  inshore  of  the  jet  center):  and  the 
shoreward  edge  of  the  surface  front  in  SST.  defined  by  the  1 1.5  C  isotherm  (blue-to-green  boundary)  during  June  and  the  12.0°C 
isotherm  during  August.  The  stars  (*)  mark  the  location  (occurrence)  of  humpback  whales  sightings  and  do  not  reflect  the  number  of 
whales  in  each  sighting. 
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Integrated  volume  backscatter  (0-100  m)  38  kHz 


Fig.  4.  The  correspondence  of  humpback  whale  sightings  (*)  with  surface  (5  m)  chlorophyll  and  volume  acoustic  backscatter  at  38  kHz 
integrated  over  O-lOOm,  during  May-June  and  July-August  2000.  A  higher  (less  negative)  value  of  integrated  backscatter  (dB) 
indicates  greater  scattering  at  that  frequency.  The  correspondence  between  humpback  whales  and  higher  chlorophyll  and  higher 
backscatter  at  38  kHz  is  evident  near  Cape  Blanco  during  spring  and  on  Heceta  Bank  during  summer. 
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the  slope  off  and  south  of  Cape  Blanco  (Fig.  6).  Of 
the  deep-diving  large  odontocetes,  Baird's  beaked 
whale  was  sighted  most  frequently  (Table  1)  and 
occurred  at  water  depths  >  1000  m  on  the  slope 
(Fig.  6).  Sperm  whales  were  only  sighted  in  slope 
waters  west  or  southwest  of  Heceta  Bank,  where 
the  bathymetry  is  steep  and  the  coastal  jet  is 
strongly  influenced  by  the  bottom  topography 
(Figs.  3  and  6).  During  August,  four  species  of 
deep-diving  large  odontocetes  (i.e.  sperm  whale, 
Baird's  beaked  whale,  Cuvier's  beaked  whale,  and 
Mesoplodon  sp.)  occurred  west  and  southwest  of 
the  submarine  bank,  in  a  band  of  surface  (5-m) 
convergence  between  the  equatorward  flowing  jet 
and  the  northward  flow  along  the  edge  of  an 
anticyclonic  mesoscale  eddy  or  meander  (Fig.  6). 
Though  velocities  were  uniformly  southward  in 
the  upper  150  m  during  June,  during  August  the 
velocity  structure  in  this  region  was  complex  with 
flow  reversals  (northward  versus  southward)  in 
both  the  vertical  and  across-shelf  structure  (see 
velocity  structure  at  Line  3,  44.24  N;  Barth  et  al., 
2005  their  Fig.  9). 

3.4.  Mesoscale  distributions  of  dolphins  and 
porpoises 


Fig.  5.  The  correspondence  of  humpback  whale  sightings  (blue 
stars)  overlaid  on  catch  distribution  for  juvenile  coho  Oncor- 
hynchus  kisutch  (red  triangles)  and  chinook  salmon  O 
tshawytscha  (red  circles)  from  surface  trawls  (Brodeur  et  al., 
2004)  conducted  during  August  2000.  Catch  data  is  overlaid  on 
surface  temperature  contours.  The  size  of  the  triangles  and 
circles  is  proportional  to  the  size  of  the  catch:  smaller  symbols 
represent  catch  of  1-5  salmon;  larger  symbols  represent  catch  of 
6-150  salmon.  Plus  signs  are  stations  sampled  where  no  salmon 
were  caught.  The  correspondence  between  humpback  whales 
and  regions  of  higher  juvenile  salmon  abundance  (near  Heceta 
Bank  and  off  Cape  Blanco)  is  clear. 


(Fig.  6).  Though  minke  whales  were  observed  on 
the  shelf,  they  were  very  rare. 

3.3.3.  Large  odontocetes 

There  is  some  suggestion  that  the  broader  slope 
at  43  43.5°N  is  less  important  to  deep-diving 
odontocetes  (or  their  prey)  than  the  steep  bathy¬ 
metry  along  a  submarine  bank,  Heceta  Bank,  and 


Pacific  white-sided  dolphin  was  the  dominant 
small  cetacean  over  the  slope  (200-2000  m)  during 
June;  however,  by  August  this  species  was  rare  in 
the  study  region  (Table  1,  Fig.  7).  Results  of  the 
multiple  logistic  regression  model  explained  44.5% 
of  the  occurrence  pattern  of  Pacific  white-sided 
dolphins  (Table  3).  A  curvilinear  response  to 
distance  to  the  shoreward  edge  of  the  jet  (geopo¬ 
tential  anomaly  2.0 m2s~2)  was  the  best  predictor 
of  the  occurrence  pattern  of  this  species.  The 
curvilinear  effect  likely  reflects  the  proximity  of 
Pacific  white-sided  dolphins  to  the  upwelling  front 
in  the  southern  region  and  their  increased  distance 
from  this  feature  in  the  northern  region  (Fig.  7). 
Near  Cape  Blanco  the  upwelling  front  and  center 
of  the  jet  extend  farther  offshore  over  the  slope, 
whereas  in  the  northern  area  these  features  are 
restricted  to  the  shelf  during  June  (Fig.  7).  Pacific 
white-sided  dolphins  occurred  typically  to  the  west 
of  the  jet  center  in  the  north,  but  were  in  the  jet  in 
the  southern  region.  The  second  most  important 
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Fig.  6.  The  distributions  of  sightings  oflarge  cetaceans  during  May-June  and  July-August  2000  relative  to  velocity  at  5  m  as  measured 
by  shipboard  Acoustic  Doppler  Current  Profiler.  A  scale  arrow  for  velocity  appears  at  right.  Symbols  for  each  species  are  as  follows: 
humpback  whale  (yellow  O);  fin  whale  (blue  O);  minke  whale  (yellow  + );  sperm  whale  (red  *);  Baird's  beaked  whale  (red  □);  Cuvier’s 
beaked  whale  (red  0);  Mesoplodon  sp.  (red  O);  and  gray  whale  (green  *).  During  August,  the  higher  density  of  deep-diving  odontocetes 
(red  symbols)  over  the  slope  in  a  region  of  surface  convergence  between  the  equatorward  jet  and  the  northward  flow  along  the  western 
edge  of  an  anticyclonic  mesoscale  eddy,  is  evident. 


variable,  and  also  the  most  important  biological 
variable  to  explain  dolphin  occurrence,  was 
integrated  acoustic  backscatter  at  38  kHz  (i.e. 
larger  targets  such  as  fish).  The  occurrence  pattern 
of  Pacific  white-sided  dolphins  was  also  signifi¬ 
cantly  positively  correlated  with  depth  and  nega¬ 
tively  associated  with  surface  salinity  (Table  2), 
reflecting  their  association  with  the  slope  and  the 
lower  salinities  occurring  there  due  to  the  influence 
of  the  Columbia  River  water  in  spring  (Fig.  7). 
Pacific  white-sided  dolphins  occurred  around  the 
edges  of  the  low-salinity  plume  in  the  northern 
region  (~44.0  44.6°N). 

Northern  right  whale  dolphins  were  often 
associated  with  schools  of  Pacific  white-sided 


dolphins  during  June.  Though  the  number  of 
sightings  was  insufficient  to  use  in  a  model  of 
occurrence  pattern,  their  distribution  on  the  slope 
was  similar  to  that  of  Pacific  white-sided  dolphins. 
This  species  was  not  observed  in  the  study  region 
during  the  July  August  survey.  Risso’s  dolphins 
(grampus)  appeared  to  be  more  widely  distributed 
and  more  abundant  in  June  than  in  August  in  the 
northern  CCS  off  Oregon  (Table  2).  They  occurred 
in  the  jet  of  higher  velocity  and  the  upwelling  front 
over  the  slope  off  Cape  Blanco  and  over  the  shelf 
and  slope  between  44.0  44.75  N  during  June. 
Though  often  associated  with  warmer  ocean 
temperatures,  this  species  occurred  in  SSTs  ran¬ 
ging  from  10-13  °C  during  June.  Their  distribution 
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during  June  and  August  suggests  a  preference  for 
the  slope  and  outer  shelf  associated  with  the 
submarine  bank,  Heceta  Bank  (44.0  44.5°N). 

The  multiple  logistic  regression  model  for  Dali's 
porpoise  explained  25.6%  of  the  variance  in  their 
distribution  during  spring,  when  they  were  rela¬ 
tively  evenly  distributed  on  the  slope  and  occa¬ 
sionally  occurred  on  the  shelf  at  depths  >100m 
(Fig.  7).  Acoustic  backscatter  at  420  kHz  (e.g., 
smaller  plankton  such  as  copepods)  was  the  most 
important  variable  to  explain  this  species'  distri¬ 
bution  during  late  May  and  early  June  (Table  3). 
There  were  no  significant  correlations  between 
Dali's  porpoise  occurrence  and  physical  variables 
during  spring.  However,  during  summer  when 
mesoscale  variability  in  the  CCS  was  higher,  the 
model  was  more  successful  and  explained  46.1% 
of  the  variance  in  occurrence  pattern  (Table  3). 
Depth  and  distance  to  the  shoreward  edge  of  the 
upwelling  front  (i.e.  geopotential  anomaly  of 
2.0  nr  s~2)  were  the  first  and  second  most  im¬ 
portant  variables  to  explain  the  variance  in 
their  distribution  during  summer.  Then,  the 
presence  of  a  warm,  low-chlorophyll,  mesoscale 
anticyclonic  eddy  or  meander  appears  to  have 
displaced  or  concentrated  Dali's  porpoises  to 
the  north  and  south  edges  of  the  meander  (Fig. 
7).  In  addition,  Dali's  porpoise  were  significantly 
negatively  correlated  with  thermocline  depth, 
reflecting  the  shoaling  of  isotherms  at  the  edges 
of  the  warm  meander.  Their  response  to  SST  was 
curvilinear,  as  Dali's  porpoise  occurred  in  a  frontal 
zone  of  intermediate  SST  (12  14CC)  between  the 
warmest  offshore  water  and  coolest  upwelled 
water  (Fig.  7).  During  summer,  Dali’s  porpoise 
occurrence  was  also  significantly  positively  asso¬ 
ciated  with  acoustic  backscatter  at  38  kHz 
(i.e.  larger  prey  such  as  fish),  which  was  the 
most  important  biological  variable  in  the  model 
(Table  3). 

Harbor  porpoise  was  the  dominant  small 
cetacean  on  the  shelf  in  both  seasons  (Table  2). 
Typically  a  coastal  denizen,  they  were  distributed 
further  offshore  (out  to  the  ^200  m  isobath)  on  the 
broad  shelf  at  Heceta  Bank  and  also  off  Cape 
Blanco  where  the  jet  and  upwelling  front  extend 
westward  (Fig.  8).  The  mesoscale  pattern  of 
harbor  porpoise  distribution  suggests  that  the 


narrow  shelf  region  between  these  features  (i.e. 
south  of  Heceta  Bank  and  north  of  Coos  Bay, 
~43.5  44.0  N)  supports  fewer  harbor  porpoises 
(Fig.  8).  The  multiple  logistic  regression  model 
explained  79.2%  and  70.1%  of  the  variance  in 
occurrence  pattern  of  harbor  porpoises  in  spring 
and  summer,  respectively  (Table  3).  Their  associa¬ 
tion  with  the  broader  shelf  or  bank  in  the  north 
makes  latitude  the  most  important  predictive 
variable  within  the  study  region  in  both  seasons. 
During  spring,  surface  salinity  was  the  second 
most  important  predictive  variable  to  explain  the 
variance  in  harbor  porpoise  occurrence.  Harbor 
porpoises  occurred  in  the  upwelled  water  of  higher 
salinity  near  the  coast.  Harbor  porpoise  occur¬ 
rence  was  also  significantly  negatively  correlated 
with  distances  to  the  shoreward  edge  of  the 
upwelling  front  (i.e.  geopotential  anomaly 
2.0  nrs-2)  and  the  center  of  the  jet  (i.e.  geopoten¬ 
tial  anomaly  2.35  m2s-2),  meaning  that  harbor 
porpoise  were  near  these  features  (Table  3  and 
Fig.  8). 

During  summer,  when  the  upwelling  front  and 
jet  occur  farther  offshore,  harbor  porpoises  are 
significantly  negatively  associated  with  only  the 
distance  to  the  shoreward  edge  of  the  surface  front 
(i.e.  12.0  °C  SST  isotherm)  and  are  significantly 
positively  associated  with  the  jet  center  (i.e., 
harbor  porpoises  are  far  from  the  center  of  the 
jet)  (Table  3  and  Fig.  8).  Harbor  porpoises  were 
associated  with  regions  of  higher  surface  chlor¬ 
ophyll  over  Heceta  Bank  and  near  Cape  Blanco  in 
both  seasons  (Fig.  8),  although  a  significant 
positive  correlation  with  value  of  the  chlorophyll 
maximum  only  occurred  during  spring  (Table  3). 
Acoustic  backscatter  at  420  kHz  (i.e.  small  zoo¬ 
plankton  such  as  copepods)  and  value  of  the 
chlorophyll  maximum  were  the  first  and  second 
most  important  biological  variables,  respectively, 
to  explain  the  variance  in  harbor  porpoise 
occurrence.  During  August,  harbor  porpoise 
occurrence  was  significantly  correlated  with  acous¬ 
tic  backscatter  at  three  frequencies  (i.e.  120,  200 
and  420  kHz);  however,  only  the  correlation  with 
120  kHz,  which  was  the  strongest  biological 
variable,  was  positive  (i.e.  harbor  porpoise  occur¬ 
rence  was  related  to  higher  backscatter  at 
120  kHz). 
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4.  Discussion 

We  were  able  to  explain  a  high  percentage  (up  to 
94%)  of  the  variation  in  cetacean  distribution 
using  measured  oceanographic  data.  Certainly, 
having  fine-scale  ocean  and  cetacean  data  on 
comparable  space  and  time  scales  was  important; 
however,  even  with  fine-scale  concurrent  sampling 
of  oceanographic  variables  and  cetacean  distribu¬ 
tion,  the  high  variance  explained  is  most  likely  due 
to  cetaceans  knowing  their  environment  and 
actively  selecting  habitat  having  certain  measur¬ 
able  and  consistent  oceanographic  qualities.  Pre¬ 
dator  knowledge  is  recognized  as  a  key  component 
in  the  relationship  between  animal  behavior  and 
ocean  structure  (Matthiopoulos,  2003).  With  our 
GLOBEC  surveys  of  the  northern  CCS,  we  were 
able  to  measure  those  important  ocean  features 
and  upwelling  structure  quantitatively  on  similar 
space  and  time  scales  to  the  cetaceans.  The  merit 
of  interdisciplinary  approaches  in  cetacean  re¬ 
search  has  been  established  through  comparisons 
of  cetacean  distributions  with  hydrographic  fea¬ 
tures  (Smith  et  al.,  1986;  Tynan,  1997;  Reilly,  1990; 
Davis  et  al.,  2002;  Jaquet  and  Gendron,  2002),  as 
well  as  associations  with  prey  (Schoenherr,  1991; 
Wishner  et  al.,  1995;  Croll  et  al.,  1998;  Fiedler  et 
al.,  1998;  Tynan,  1998;  Benson  et  al.,  2002;  Gill, 
2002).  However,  the  necessity  in  some  studies  to 
pool  data  into  large  ^2°  latitude-longitude  blocks 
can  impose  a  minimum  level  of  resolution  in 
models  (Jaquet  et  al.,  1996;  Reilly,  1990).  Com¬ 
parisons  of  mammal  distributions  with  monthly 
mean  indices  of  upwelling  (Keiper  et  al.,  2005)  also 
can  negate  the  importance  of  upwelling  dynamics 
on  finer  temporal  scales.  In  the  present  study,  the 
synoptic,  high  resolution  of  the  biophysical 
processes  associated  with  upwelling  and  the 


circulation  of  the  northern  CCS  provided  a  close 
spatial  and  temporal  comparison  with  the  cetacean 
survey  data. 

Flow-topography  interactions  between  the 
northern  California  Current  and  Heceta  Bank 
(ca.  44.0  44.5CN)  strongly  influenced  cetacean 
distributions.  The  ADCP-derived  current  vectors 
suggest  a  recirculation  over  the  bank,  with  south¬ 
ward  flow  over  the  northern  and  western  regions 
of  the  bank  and  some  eastward  and  northward 
flow  at  the  southern  edge.  Humpback  whales 
occurred  on  the  bank  during  August,  when  their 
distribution  was  significantly  correlated  with  high 
chlorophyll  concentration  and  higher  backscatter 
at  38  kHz  (i.e.  larger  prey  such  as  fish).  The  highest 
catch  of  juvenile  chinook  Oncorhynchus  tsha- 
wytscha  and  coho  salmon  O.  kisulch  (10-25  fish 
per  catch)  (Brodeur  et  al.,  2004),  as  well  as  the 
highest  abundance  of  Pacific  sardine  Sardinops 
sagax  (1000  10,000  fish  per  catch;  R.D.  Brodeur, 
unpublished  data  NMFS)  occurred  where  hump¬ 
backs  concentrated  on  Heceta  Bank.  During 
August,  the  highest  densities  and  largest  daytime 
patches  of  euphausiids  ^10mm  in  length  (domi¬ 
nated  by  E.  pacified)  also  occurred  over  Heceta 
Bank  (Ressler  et  al.,  2005). 

Therefore,  it  appears  that  by  August,  humpback 
whales  were  responding  to  a  cascade  of  trophic 
dynamics  enhanced  by  the  flow-topography  inter¬ 
actions  and  the  broad  upwelling  signature  over  the 
bank.  The  extraordinary  success  of  our  multiple 
logistic  regression  model  to  explain  variance  in 
humpback  occurrence  during  August  (94.4%),  is 
likely  due  to  the  success  of  the  models  to  capture 
the  occurrence  of  humpbacks  near  these  ‘hot 
spots'  of  increased  prey  density. 

Although  Heceta  Bank  is  among  the  largest 
submarine  banks  on  the  west  coast  (Hickey,  1998), 
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Fig.  7.  The  correspondence  of  sightings  of  Pacific  white-sided  dolphins  (O)  and  Dali's  porpoises  (*)  with  sea  surface  (5  m)  temperature, 
surface  salinity,  and  geopotential  anomaly  at  5  m  relative  to  100  m,  during  May-June  and  July-August  2000.  Several  features  related  to 
upwelling  were  included  in  the  models  of  cetacean  occurrence:  the  center  of  the  equatorward  jet,  defined  by  geopotential  anomaly 
2.35  nrs  ~  (the  double  heavy  line);  the  shoreward  edge  of  the  jet,  defined  by  geopotential  anomaly  2.0  m2s-2  (first  single  heavy  line 
inshore  of  the  jet  center);  and  the  shoreward  edge  of  the  surface  front  in  SST,  defined  by  the  1 1.5°C  isotherm  (blue-to-green  boundary) 
during  June,  and  the  12.0  C  isotherm  during  August.  The  symbols  mark  the  location  of  sightings  and  do  not  reflect  the  number  of 
dolphins  or  porpoises  in  each  sighting.  The  occurrence  of  Pacific  white-sided  dolphins  in  or  west  of  the  center  of  the  jet,  in  warmer, 
lower  salinity  surface  water,  is  evident  during  June.  During  August,  the  occurrence  of  Dali's  porpoise  between  the  center  and 
shoreward  edge  of  the  jet  and  upwelling  front  is  clear. 
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internal  hydraulic  flows  associated  with  even 
smaller  features  on  this  bank  (i.e.  Stonewall  Bank) 
can  greatly  enhance  turbulence  in  the  bottom 
boundary  layer  (Nash  and  Mourn,  2001).  Fine- 
scale  process  studies  to  examine  associations 
between  foraging  cetaceans,  their  prey  and  the 
circulation  of  banks  warrants  further  study.  In  the 
southeast  Bering  Sea,  foraging  whales  occur  on  the 
middle  shelf  (50  100  m)  where  the  mean  currents 
are  sluggish,  but  tidal  currents  are  energetic 
(Tynan,  2004).  Here,  we  hypothesize  that  the 
enhanced  vertical  and  horizontal  mixing  asso¬ 
ciated  with  Heceta  Bank  is  linked  to  higher  prey 
densities  and  hence,  improved  foraging  for  ceta¬ 
ceans,  especially  humpback  whales  and  harbor 
porpoises.  During  August,  humpback  whales 
were  observed  feeding  along  and  near  a  series  of 
parallel  alongshore  fronts  on  Heceta  Bank.  The 
fronts  may  have  resulted  from  internal  wave  field 
passage. 

The  surveys  during  June  and  August  2000 
provided  an  opportunity  to  compare  the  seasonal 
shift  in  the  position  of  the  equatorward  coastal  jet 
with  the  distribution  of  cetaceans.  Typically,  the 
jet  forms  near  the  coast  during  spring  upwelling 
and  moves  offshore  with  sustained  upwelling 
(Mourn  and  Nash,  2000).  The  longitudinal  posi¬ 
tion  of  the  shoreward  edge  of  the  surface  front 
however  is  not  greatly  different  between  June  and 
August  2000,  although  during  August  the  meso- 
scale  structure  was  enhanced  and  the  surface  front 
was  especially  convoluted  on  Heceta  Bank.  In 
addition,  the  intrusion  of  a  warm  low-chlorophyll 
mesoscale  meander  offshore  (ca.  43  44  N)  during 
August  restricted  the  offshore  movement  of 
upwelled  waters  and  compressed  the  system  to 
the  north,  south  and  east.  This  feature  also  may 
have  compressed  distributions  of  slope-inhabiting 
cetaceans,  such  as  Dali's  porpoise,  to  the  periphery 


of  the  meander  between  the  center  and  shoreward 
edge  of  the  equatorward  jet.  The  increased 
mesoscale  variability  during  summer  in  this  way 
improved  the  percentage  of  explained  variance  in 
Dali's  porpoise  occurrence  pattern.  Shoaling  iso- 
pycnals  at  the  edge  of  the  meander  may  also  have 
compressed  the  depth  range  of  their  fish  or 
cephalopod  prey  (Walker,  1996),  thereby  enhan¬ 
cing  foraging  at  the  periphery. 

Pacific  white-sided  dolphin,  which  was  the 
dominant  small  cetacean  over  the  slope  in  June, 
was  very  rare  in  the  system  during  August; 
however,  their  high  seasonality  likely  reflects  a 
movement  or  migration  from  the  region  more  than 
a  regional  displacement  by  the  warm  mesoscale 
meander  in  August.  The  strong  correlation  be¬ 
tween  Pacific  white-sided  dolphins  and  backscatter 
at  38  kHz  suggests  that  during  spring  they  were 
following  mesopelagic  fish,  although  cephalopods 
are  also  important  in  the  diet  (Walker  and  Jones, 
1990). 

The  distributions  of  harbor  porpoises  and 
humpback  whales  were  significantly  associated 
with  the  position  of  the  alongshore  upwelling 
front  in  both  seasons;  however  more  humpback 
whales  occurred  offshore  on  the  western  edge  of 
the  jet  and  upwelling  front  in  June.  In  contrast,  by 
August,  humpback  whales  and  harbor  porpoises 
concentrated  in  high-chlorophyll  shelf  waters  on 
the  inshore  side  of  the  jet  and  surface  front  at 
Heceta  Bank  and  near  Cape  Blanco.  The  region 
inshore  of  the  surface  front  is  associated  with 
convergence  (i.e.  cross-shelf  streamlines  leave  the 
surface  shoreward  of  the  density  front;  Allen  et  al., 
1995,  their  Figs.  4  and  5),  and  therefore,  with  the 
potential  concentration  of  prey.  The  flow-topo¬ 
graphy  interactions  that  contribute  to  the  posi¬ 
tions  of  the  jet  and  the  inshore  edge  of  the 
upwelling  front  are  here  recognized  as  key 


<- 

Fig.  8.  The  correspondence  of  sightings  of  harbor  porpoises  (*)  with  surface  (5  m)  chlorophyll,  surface  salinity,  and  geopotential 
anomaly  at  5  m  relative  to  100  m,  during  May-June  and  July-August  2000.  Several  features  related  to  upwelling  were  included  in  the 
models  of  cetacean  occurrence:  the  center  of  the  equatorward  jet,  defined  by  geopotential  anomaly  2.35  m2s~2  (the  double  heavy  line); 
the  shoreward  edge  of  the  jet,  defined  by  geopotential  anomaly  2.0  m2s-2  (first  single  heavy  line  inshore  of  the  jet  center);  and  the 
shoreward  edge  of  the  surface  front  in  SST,  defined  by  the  1 1.5  °C  isotherm  (blue-to-green  boundary)  during  June  and  the  12.0  °C 
isotherm  during  August.  During  June,  harbor  porpoise  occurrence  is  restricted  to  the  shelf,  near  the  upwelling  front  and  shoreward  of 
the  center  of  the  jet.  The  high  density  of  harbor  porpoise  in  upwelled  waters  on  Heceta  Bank  during  August  corresponds  with  the 
shoreward  edge  of  the  surface  front,  high  surface  chlorophyll,  and  high  salinity. 


164 


C.T.  Tynan  et  al.  /  Deep-Sea  Research  II 52  (2005)  145-167 


processes  to  include  in  predictive  models  of 
cetacean  occurrence  in  the  CCS.  Similar  processes 
and  features  may  also  be  important  in  determining 
cetacean  distributions  in  other  boundary  current 
systems  of  the  world. 

Humpback  whales  arrive  in  the  northern  CCS  in 
spring  after  spending  the  winter  off  Mexico 
(Calambokidis  et  al.,  2000).  During  spring,  hump¬ 
back  whale  occurrence  was  significantly  correlated 
with  the  warmer  SSTs  over  the  slope;  however, 
since  some  humpbacks  occurred  in  the  cold 
recently  upwelled  water  of  higher  backscatter  at 
38  kHz  (i.e.  likely  fish),  it  is  assumed  that  hump¬ 
backs  can  tolerate  the  cold  (~9°C)  conditions 
after  they  arrive.  Humpbacks  did  appear  to  avoid 
the  lowest  salinity  water  of  Columbia  River  origin; 
however,  they  occurred  at  the  salinity  front 
associated  with  the  offshore  plume.  During  spring, 
humpbacks  off  and  south  of  Cape  Blanco  were 
significantly  correlated  with  the  center  of  the 
equatorward  jet.  Our  models  demonstrate  the 
important  predictive  capacity  of  dynamic  height; 
it  is  encouraging  that  this  variable  could  be  applied 
in  models  of  occurrence  in  other  boundary  current 
systems. 

Both  slope-associated  Dali's  porpoise  and  shelf- 
associated  harbor  porpoise  were  significantly 
associated  with  integrated  backscatter  at  420  kHz 
during  spring.  For  both  species,  backscatter  at 
420  kHz  was  the  most  important  predictive  biolo¬ 
gical  variable  during  June.  It  is  possible  that 
backscatter  of  zooplankton  at  the  high  frequencies 
serves  as  a  better  proxy  for  the  porpoises'  prey  (i.e. 
fish  and  cephalopods),  especially  for  species  that 
may  diurnally  migrate  from  depths  below  100  m 
(or  deeper  than  the  depth  integration  of  acoustic 
backscatter,  0-1 00  m).  However,  though  it  is 
tempting  to  assume  that  backscatter  at  420  kHz 
is  indicative  of  copepods,  an  important  caveat  in 
the  interpretation  of  integrated  backscatter  is  that 
high  abundances  of  gelatinous  zooplankton  can 
contribute  significantly  to  backscatter  measured  at 
^200  kHz  (Monger  et  al.,  1998).  Therefore, 
interpretation  of  the  associations  between  ceta¬ 
ceans  and  backscatter  at  higher  frequencies  (i.e. 
200  or  420  kHz)  awaits  further  resolution  of  the 
bioacoustic  data  to  more  specific  taxonomic 
definition  using  net  collections. 


The  presence  of  sperm  whales  in  slope  waters  off 
Heceta  Bank  in  spring  and  summer  suggests  that 
their  cephalopod  or  fish  prey  are  predictably 
abundant  over  the  steep  and  rough  topography 
in  this  region.  In  general,  large  cephalopods  and 
fish  are  the  primary  and  secondary  prey,  respec¬ 
tively,  of  sperm  whales  off  California  (Kawakami, 
1980).  There,  the  large  squid  Moroteuthis  robusta 
is  an  important  prey  species  of  sperm  whales 
(Fiscus  et  al.,  1989);  specimens  of  >  1  m  lengths 
have  been  caught  off  Oregon  in  water  deeper  than 
150  m  (Pearcy,  1965).  Though  squid  were  not 
specifically  studied  during  GLOBEC,  some  earlier 
studies  suggest  that  an  influx  or  migration  of  squid 
into  slope  waters  off  Oregon  occurs  during 
summer  (Pearcy,  1965). 

The  association  of  sperm  whales  and  beaked 
whales  (i.e.  B.  bairdii ,  Z.  cavirostris ,  and  Mesoplo - 
don  sp.)  with  a  dynamic  region  of  surface 
convergence  between  the  southward  flowing  jet 
and  the  northward  flow  along  the  edge  of  an 
anticyclonic  mesoscale  eddy  during  August  sug¬ 
gests  that  flow-topography  interactions  along  the 
slope  near  Heceta  Bank  are  important  to  large 
odontocetes  and  their  prey.  We  hypothesize  that 
strong  turbulence  associated  with  rough  topogra¬ 
phy  (Kunze  and  Toole,  1997)  and  interactions  of 
the  internal  tide  with  a  critical  slope  (Thorpe  et  al., 
1990;  Mourn  et  al.,  2002)  may  influence  foraging 
sperm  whales  and  beaked  whales  through  linkages 
between  their  cephalopod  or  fish  prey  and 
enhanced  nepheloid  layers  at  regions  of  high 
mixing  over  the  continental  slope. 
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